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A novel quantum interference effect in disordered quasi-one-dimensional rings in the inhomoge- 
neous magnetic field is reported. We calculate the canonical disorder averaged persistent current 
using the diagrammatic perturbation theory. It is shown that within the adiabatic regime the av- 
erage current oscillates as a function of the geometric flux which is related to the Berry phase and 
the period becomes half the value of the case of a single one-dimensional ring. We also discuss the 
magnetic dephasing effect on the averaged current. 

PACS numbers: 03.65.Bz, 73.23.Ra, 73.10.-d 

Since the discovery of the Berry phase, |jj there has been much interest in the study of topological effects in the 
fields of quantum mechanics and condensed matter physics. Q The typical example to illustrate the Berry phase is 
the Aharonov-Bohm (AB) effect || in the mesoscopic ring, where the relative phase would accumulate on the wave 
function of a charged particle due to the presence of a electromagnetic gauge potential. Similarly, when a quantum 
spin follows adiabatically a magnetic field that rotates slowly in time, the spin wave function acquires an additional 
geometric phase (Berry phase) besides the usual electromagnetic phase in the static magnetic fields. 

Loss et al. [Q studied the persistent current in an ideal one-dimensional ring in the presence of a static inhomo- 
geneous magnetic field by examining the coupling between spin and orbital motion through the Zeeman interaction. 
Using the imaginary time path integral method, they showed that the spin wave function accumulates the Berry phase 
when the spin of an electron traversing an AB ring adiabatically follows a textured inhomogeneous magnetic filed 
with a tilt angle and this phase leads to persistent equilibrium charge current. 

In this paper we investigate the persistent charge current of the quasi-one-dimensional disordered rings Q (see the 
inset in Fig. 2) and compare it with the results of Loss et al. Q We compute the disorder averaged persistent current 
with the help of the diagrammatic technique and find that the average current is given by the sum of the term which 
oscillates as a function of the Berry phase and the term which does not depend on the the Berry phase. We show that 
the latter term is negligibly small compared with the former term in the adiabatic regime, so that we will be able to 
clearly observe the Berry phase dependent persistent current experimentally. 

We begin by considering a quasi-onc-dimensional ring of circumference L x = 2irr and volume V = L x L y L z . The 
ring is embedded in an static inhomogeneous magnetic field B. For a spin-1/2 electrons of mass m and charge e, the 
system may be described by the Hamiltonian 
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-A e ™(r) + U (r)-- 3MB B(r).<T, (1) 



where p, r, g, /ib, and ha/2 are the momentum, position, g factor, Bohr magneton, and spin, respectively. The 
operator u(r) represents the spin-independent random impurity potential and, A em is the electromagnetic gauge 
potential, with B = V x A em relating it to the magnetic fields. In the following we specialize in the case of 
inhomogeneous magnetic fields with constant magnitude B, and we have parametrized B in terms of the spherical 
polar angles \ and r\ so that it has Cartesian components B (sinx(r) cos rj(r), sinx(r) sin ^(r), cos x(r)), with the angles 
X and rj being smooth functions of position. 

Using the Green's function, the thermal average of the persistent current is given by 

/= / ^"E E /x/(e)[< Q (k,k, E )-G« a (k,k, e )], (2) 

J k Q=±l 

where a = ±1 is the spin index and /(e) is the Fermi-Dirac distribution function. The current vertex is given by 
I x = —ehk x /mL, where k x is the x component of the wave number vector. Moreover is the retarded (advanced) 

Green's function and it denotes the exact Green's function before impurity averaging. To determine the canonical 
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disorder-averaged persistent current, we follow the procedure described in Refs. [6~8] The chemical potential is set 
to /i = no + 5p and we expand Eq. (||) to first order in Sp. Then we obtain the averaged persistent current at a fixed 
particle number as follows: 

( 7 ( $em )) " "f 9§^({ E^O*, *""*)} ) , (3) 



where $ em is the electromagnetic flux and A is the mean level spacing. The fluctuation of the particle number 
{(12a <^«) 2 ) i s related to the two-point correlator of the density of states p a , |J 

(I E SN a (fx ,^ em )\ \ = V 2 r dex ^ de 2 f(e 1 )f(e 2 )J2K a ^(e 1: e 2 ), (4) 

\ U=±l J / J -°° J -°° a,a< 

where 

K a>a '(ei,E2) = (p a {ei)p a '(S2)) - (Pa(£l)) (Pa'(£2)} 

-^F2 / dr i / rfr 2 Re[(G^ a (r 1 ,r 1 ;£ 1 )G^ )a ,(r 2 ,r 2 ;e 2 )>-(G^ Q (r 1) r 1 ;e 1 )>(G^ )a ,(r 2) r 2 ;e 2 )>]. (5) 

In order to calculate K a ^ a i , we use the diagrammatic perturbation method JlO| which is a powerful tool in the theory 
of disordered system. In the language of diagrammatics, the dominant contribution to K a%a * arises from the exchange 
of the two Cooperon (particle-particle) ladders between two closed loops |||| (see Fig. 1). On the other hand, the 
Diffuson process gives no contribution to the average current because it has no dependence of the magnetic flux 
$™. Thus we need only the Cooperon process in order to calculate the average persistent current. As usual, we are 
now interested only in the field dependence resulting from interfering paths. Therefore we ignore the less important 
field dependence of the averaged single Green's function. The two-point correlator then translates into the following 
expression: 

K a , a >{e 1: e 2 ) = — -jRe / dxi / dx 2 C a>a > (xi,x 2 ;£i - ei)C a ^ (x 2 ,x 1 ;ei - e 2 ) . (6) 
l-K z V z n J J 

In this equation we have used the definition of the particle-particle pair propagator 

Ca t a> (xi,X 2 ;£l - £ 2 ) = (G£ Q (x2,Xi;£i)G^ Q /(x2,Xi;£2)) , (7) 

where p(0) is the density of states (per unit volume and spin) at the Fermi surface. Following Ref. [11] we have 
evaluated the pair propagator C a a i by using the quasi-classical Feynman path-integral method JT3] and obtained 
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FIG. 1. An example of diagrams showing the leading contribution to the average persistent current. 
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C Q , Q /(xi,x 2 ;ti,t 2 ) =e(t 2 -h) / CRexp(-i5^[R]), (8) 

jR(t 2 )=x 2 V 7 

where the effective action is given by 
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[R] = - 4^ / dtn'-^-l dtA- A em (R(t)) - ^ (i a - *!)(a - a') - 7a [R] - 7<*' [R] ■ (9) 
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Here the Berry phase 7 a [R] which arises from adiabatic approximation for the spin (dynamical Zeeman) propagator 
is of the form 

7q [R] = f 2 dtR.A£(R), (10) 

where A^(R) is the spin dependent geometric gauge potential and is given by A^(R) = (a/2)Wi] [cos\ — 1]- It should 
be note that the expression Eq. (@) is only valid in the adiabatic regime in which the spin of the electron adiabatically 
follows the local direction of the non-uniform magnetic field. This adiabaticity requires that the precession frequency 
lob = g[J,BB/2h is large compared to the reciprocal of the diffusion time = L x /D (D is the diffusion constant) 
around the ring, i.e., ujsTd ^> 1, or equivalcntly B » B c = 2Eth/ 9^B- Pi - pit 

The path-integral representation for the pair propagator can be transformed into a differential equation, 



7 + D j-«T^7 - A « ( x - A «' ( x ')| - iu B (a, ~ a') 



C QjCe < (xi,x 2 ;ti,t 2 ) =S(xi- x 2 ) S (ti - t 2 ) , (11) 



which is characterized by the presence of the spin-dependent gauge potential A a = (e/h)A e,n + A^. We now evaluated 
the pair propagator for the case of the metallic rings, having height L z and thickness L y . The magnetic field was 
assumed to have constant magnitude, and not to vary appreciably either radially across the wall or parallel to the 
ring axis. Thus the gauge potential A„ can be replaced by its tangential component (<j) is the azimuthal angle). 
This can be expressed as 2-ktA^ = Q a = <& em /<I>o + a$ 9 , where $ em is the electromagnetic flux through the area nr 2 
and the geometric flux is given by 

$fl = ^/ «ty[cosxW)-l]0*»?(0). (12) 

The geometric flux take a simple form in the case of a cylindrically symmetric texture (rj = <fi), i.e., <& 9 = (cos^ — l)/2, 
where \ is the constant tilt angle away from the cylinder axis. According to the standard method, [fl4[ we can show 
that the pair propagator C a ,a' is explicitly given by [pT[ 



1 v-^ f°° duj 

C a .a' (Xi, X 2 ; ti, t 2 ) — — > I —C a . a > (k, t 
V J-ao 2TT 



w \ e ik-(xi-x 2 )-iu)(ti-t2) ) Mg\ 

C a ^ (k, w) = j , (14) 

- m a , a >) + — + Dk\ + ^[n x - ($ Q + $ Q -)] 2 

where k = (27m;,; kj_), k^ = (2Ttn y /L y ,2Tm z /L z ) (n Xl n y ,n z are integer) and w Q Q ,< = ws(a-a'). For convenience 
we have incorporated the effect of dephasing 1/t v = D/L 2 ^ (L v is the phase coherence length) directly into the pair 
propagator. Note that for diagonal component of the pair propagator C Q . a , the Zeeman terms cancel, i.e., w a ^ a = 0. 

In the following we shall calculate the diagonal part (a' — a) and the off-diagonal part (a' — —a) of the average 
current separately, namely, 

(/($ em )) = (l D ($ em )) + (l OD (<p em )) . (15) 

The diagonal part is attributed to the interference between same spin states. If the thickness and the height of the 
ring is small compared with the length L v , then the integration over kj^ should be replaced by a summation with only 
the term corresponding to k^ = retained. Substituting the Fourier transform of Eq. ( p^ ) into Eq. (^|) and using the 
Matsubara method, therefore, we obtain the diagonal part of the average current 



3 



a=±l 



" " " " v 1 vi + — + 4ir 2 E Th (n x + 2$ Q ) 2 



(16) 



where j3 — 1/kgT and E-ph = hD/L 2 is the Thouless energy and i/g = 2tt£/(3 (£ is an integer) is the boson Matsubara 
frequency, respectively. In the limit of zero temperature, the Matsubara sum turns into an integral — ► 2tt / ' (3 J dv, 
which is easily evaluated. This yields for the diagonal part of the averaged persistent current at T — as 



(l D {* em )) = Tt E Eexpf-n^)sin(4^ci> a ) 



a— ±1 n— 1 



(17) 



Iq = evp/L x is the current carried by a single electron state in an ideal one-dimensional ring and M — k 2 F V/ L x is the 
effective channel number, where vp is the Fermi velocity and kp = mvp/fi is the Fermi wave number. In Eq. (FT 



the magnetic dephasing length L B is given by l/(L B ) 2 = 1/{L V ) 2 + \/(L B ) 2 and L B = v / 3 < I > o/(27ri a \B Z \) accounts 
for the dephasing effect of the magnetic field penetrating into the sample region, Jl5[ where B z is the z component of 
B. Therefore we will observe the persistent current oscillation with somewhat reduced amplitude. |L6| It is important 
to note that unlike the case of the one dimensional ring Q] the inhomogeneous magnetic fields leads to the dephasing 
effect to the persistent current. 

For the off-diagonal part of the persistent current which results from the interference between opposite spin states, 
same procedure as above yields 



(/ 0£l (<j> em )) 
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(18) 
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(20) 



where T = D /[2u>b(L b ) 2 ] and £ = y/hujE/Erh- Note that the off-diagonal term does not depend on the Berry 
phase. In contrast to the diagonal term, we can expect that this term gives negligibly small contributions to the 
average current in the adiabatic regime because the Cooperon pole has a non-zero Zeeman term. Figure 2 shows 
the B dependence of (l D ) and (l OD ) of quasi-one-dimensional rings (L x = 11/im and L y = 2.0/im) in the case of 
the cylindrically symmetric texture. As expected above, (J 0£) ) is much smaller than (Z 13 ) in the adiabatic regime, 
B ^> B c ps lOGauss. Therefore the average persistent current is given by 



</($ em )) « (/ D ($ em )) = A ^ 



sin (47r$ Q ) 



Q =± 1 cosh 



(21) 



cos (47r<i> Q ) 



which is a periodic function of & em and $ 9 with period $0/2 and 1/2, respectively. These periods are half the value 
of the case of a single one-dimensional ring. [Q As the half flux periodicity of the electromagnetic flux in disordered 
rings, Ji| that of geometric flux is ascribed to the ensemble averaging for fixed particle number: averaging eliminates 
the first Fourier component of the current although the second components which results from the interference between 
time-reversed trajectories survives. If the magnetic field is homogeneous then the Berry phase vanishes, i.e., $ 9 = 0. 
In this case the expression Eq. (|2l|) is reduced to the well-known result. ||p| 

In summary, we have investigated the persistent current of disordered rings in the inhomogeneous magnetic field. 
Halving of the geometric flux period with respect to the single ballistic one-dimensional ring H is obtained in the 
canonical disorder-ensemble average. We have also showed that the existence of magnetic field which penetrates the 
sample causes the dephasing effect. 

Finally we point out that the phenomena investigated in this paper will also appear in doubly-connected chaotic 
billiards, pqj We expect that these systems will open a new area to explore the effect of chaos on the Berry phase. 
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FIG. 2. Magnetic field dependence of the persistent current of quasi-one-dimensional disordered rings embedded in the 
cylindrically symmetric texture with a fixed value of $ em /$o = 1/8 and \ — n/3. Plotted is the result for D = 9.0 x 10~ 3 m 2 /s 
and L v = 5.5^m. [17] The solid (dotted) line represents the diagonal (off-diagonal) part of the persistent current. Inset shows 
a schematic drawing of a disordered ring in inhomogeneous magnetic fields (arrows). 
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